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We generalize the hybrid magneto-fluid model of a charged fluid interacting with an electromag- 
netic field to the dynamics of a relativistic hot fluid interacting with a non-Abelian field. The fiuid 
itself is endowed with a non-Abelian charge and the consequences of this generalization are worked 
out. Applications of this formalism to the Quark Gluon Plasma are suggested. 
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I. INTRODUCTION 



Recent experiments at the relativistic heavy ion col- 
lider (RHIC) have shed light on the behavior of hot 
dense nuclear matter 111. The conjecture, that quarks 
and gluons de-confine and become a plasma in extreme 
conditions 0, is close to being experimentally proved. 
However, it has been realized that the de-confincd quark 
gluon matter that has been revealed at RHIC, far from 
being the weakly interacting coUisionless plasma envi- 
sioned by theorists, is, in fact, behaving more like a quark 
gluon liquid, or a strongly interacting plasma 0. The 
QGP liquid (or strongly interacting plasma) is dense, but 
seems to flow with very little viscosity. It flows so freely 
that it approximates an ideal, or perfect fluid, the kind 
governed by the standard laws of hydrodynamics. Thus, 
most of the phenomenological input for the explanation 
of the data at RHIC comes from adopting a hydrody- 
namic approach to the plasma and consequently, fluid dy- 
namics is used to deduce its properties. For a proper fluid 
dynamical description of the quark gluon fluid the fact 
that quarks and gluons have a non-Abelian charge has 
to be taken into consideration. There have been theories 
focusing on various aspects of an ideal fluid in interaction 
with Yang-Mills fields called Yang Mills rnagnetohydro- 
dynamics or chromohydrodynamics 0,|1,IEIMM0| ■ The 
relativistic heavy ion experiments coupled with recent 
studies by Jackiw and his collaborators on the "Clebsch 
representation" of Yang Mills fluid dynamics have lead 
to a resurgence of interest in these theories P, ^3 ■ 
In tune with this revival, we investigate the dynamics 
of a relativistic hot fluid with a non-abelian charge in 
terms of a model which unifies the Yang Mills field with 



the flow field strength tensor j^. Apart from its pos- 
sible phenomenological applications, the motivation for 
the following work is based on the aesthetic criterion of 
unifying the fluid field and the Yang Mills field into a 
Yang Mills " magneto- fluid" by a " gauge principle" . The 
similarities between one gauge field theory - electromag- 
netism and fluid dynamics have been explored extensively 
and fluid flow has been shown to have a formal equiva- 
lence with a gauge theory [Hill. It has also been shown 
that a fully antisymmetric flow tensor, resembling the 
electromagnetic field, can be constructed and the unifi- 
cation is achieved by defining an effective field strength 
tensor that combines app ropriately weighted electromag- 
netic and flow fields [T3|. Is a consistent and useful non- 
Abelian generalization of this genre of flow-field unifica- 
tion possible? This investigation constitutes the theme 
of this paper. 



II. ABELIAN MAGNETO-FLUID 
UNIFICATION. 



First, let us recapitulate the salient features of 
Abelian(Maxwell) Magneto-fluid unification j^. Al- 
though Maxwells electrodynamics provides equations of 
motion for the electric and magnetic fields, for describing 
their interaction with matter fields (charged particles), 
the Lorentz force law has to be independently postulated. 
In contrast, in a gravity coupled plasma, a natural conse- 
quence of general covariance is the conservation of energy 
and momentum, and the Lorentz force law for charged 
particles moving in a gravitational field |0| can be de- 
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rived from the field equations to be : 

VuU = -^U.F, (1) 
m 

where, U is the tangent vector to a geodesic (the velocity 
vector). In the limit of a weak gravitational field (the 
flat space-time limit, — > 9^ ), the component form of 
Eq.(l) reads : 

WdM'' = —WF^ (2) 
m 

Because of the antisymmetry of F^^, , contraction with 
reduces the right hand side of Eq.(2) to zero. Thus, with 
no loss of generality, the coefficient of If^, on the left hand 
side of the equation, should also be anti-symmetrized. 
Since C/^, the four- velocity, obeys U^U^ = —1 (implying 
U^'^yU^' = 0), Eq.(2) may be written as 

U'^dM'' - V^d^Un, = —U^'F\ (3) 

TO 

With the definition 

P^u =d^U,- d,U^. (4) 
we can write Eq.(3) as 

77? 

t/"(-Pp. + ^^M-) = 0. (5) 

For deriving equations of motion for point particles, usu- 
ally a limiting procedure is invoked 15]. For the motion 
of fluids, however, a small volume element of the fluid is 
the limiting element and the statistical properties of the 
fluid come into play. Recently suggested that the 
Sfii, = d^fUv — dvfU^ must replace the particle P^^, for 
a new and natural "minimal coupling" to describe a fluid 
interacting with the Maxwell fields, where / represents a 
temperature dependent statistical attribute of the fluid, 
and is related to the enthalpy h, number density n and 
mass TO of the fluid by the relation h = mnf{T). In terms 
of S'^y, the "fluid Lorentz equation" derived in is 

m 

Td^a = qiF^"" + -S^'')U^ (6) 

where, a is an entropy density. The limiting procedure 
to the point particle case, then, is simply equivalent to 
/C/M ^ u^'^ S^, ^ as / ^ 1 (T ^ 0). 

The curvature F^i, corresponding to the connection 
is obtained from the commutator of two covariant deriva- 
tives = dfj, — iqA^. In a similar vein, we can define 
a "unified" connection Q'^ = + ^fU'^, which corre- 
sponds to a minimally coupled hot magnetofluid. The 
new Abelian covariant derivative for the unified field 

Df, =df_,- iqQf^ ^ df_, ~ iqA^ - imfU^. (7) 

leads to the unified curvature 

[Df^, Dy] = -iqF^i, - imS^„, (8) 

where F^^, = df_,A^ - d^A^ and 5*^^ = d^^fU^, - d^JU^ 

This new, temperature dependent minimal coupling 
procedure has many interesting consequences which are 
explored in 0,0. 



III. THE YANG MILLS MAGNETO-FLUID 
TENSOR 

Now we return to the main theme of this work-the 
dynamics of a non- Abelian fluid. The non- Abelian gauge 
field is represented by 

F\u - d^A% - d,A\ - ig[A^, A,]\ (9) 

where, = iC^^A^ ^.A"^, Cbc are the structure 

constants of the gauge group and g is the gauge " charge" . 

Since F now has a gauge index, the right hand side 
of the Abelian equations of motion suggests a general- 
ization of the fluid flow vector to include a gauge, or a 
non-Abelian index. The RHS of the equations Q can 
then be written as : {q/m)U^aF'^'^ f_i. Correspondingly, 
the LHS, U^df^U''' requires a non-Abelian generalization. 
This mandates giving the flow field C/^ a non-Abelian in- 
dex, and we are led to a generalization of the Abelian 
flow tensor 5*^,^ to S'^^^. 

Following the Abelian route, an explicit form of S""^^ 
is given by evaluating the non-Abelian curvature for the 
generalized non-Abelian covariant derivative, 

D^ = d^-ig[A^,]-im[fU^,]. (10) 

The corresponding curvature 

[D^,D,r ^~igF^,-zmS-^, (11) 

coupled with Eq.9, defines the YM generalization of 

S% = dJU\-d,fU%-tmf[U^,U,r 

- igf[A^,U,r^igf[U^,A,r. (12) 

written, succintly, as 

= T^'^im) - T^'im) ~ ^mf[u^, un (la) 

where = 9^ — ig[A^, ] is the ordinary non-Abelian 
gauge covariant derivative. Notice that Sj^" encompasses 
the pure flow-field as well as the interaction. 

In the presence of a matter gauge current, the Yang 
Mills Field evolves as 

The L.H.S of Eq.(14) is easily related to the energy 
momentum tensor of the gauge field 

through the Bianchi identity 

q^upap ^ -j;)app^ _^ -pP p'^'^ ^ Q, (16) 

leading to 

d'^e.p = -F^^V.FP^ = (17) 
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The evolution of Yang Mills potential Af^ = ^I^Ta is dic- 
tated by the matter current in addition to the inherent 
nonlinearities in field equations. 

In the conventional particle description, the current 
is constructed from the wave function, which itself 
evolves covariantly in the gauge field. Ideally, for strongly 
interacting matter such as a quark gluon plasma, such a 
current has to be constructed from a collective many- 
body wave function. This procedure is very cumbersome 
and sometimes not very illuminating. Thus, a fluid de- 
scription, in terms of flow and thermodynamic variables, 
is useful, because, it captures a complicated dynamics 
in terms of a few collective variables. We are "viewing" 
a strongly interacting many particle system through a 
set of representative "flow" fields. We consider not just 
one but several flow flelds (labeled by a species index s) 
denoting the particles (quarks, anti-quarks, gluons) in- 
teracting with the non-Abelian gauge field. Each species 
can in principle have different charges, densities, temper- 
atures etc. Each species is taken to be a perfect fluid 
with an energy momentum tensor of the form 

rr = PsV''" + hsU^^^U"'''', (18) 

where, Ps is the pressure, and the enthalpy density hs is 
given by 

K^{ms/T) 

hs = msnR^s—— — = msnR^sJs[T), (19) 

K2{ms/T) 

and n/j.s measures the density in the rest frame for the 
given species. 

Each of the non-Abelian species contributes a flux 

r^a,. = rifl^C/,^,,, (20) 

towards the total non-Abelian current 

•^a=E5«rV., (21) 

s 

where gs is the gauge coupling for the species 's" . 
The matter fields contributing to evolve covariantly 
under the action of the non-Abelian covariant deriva- 
tive. The multi-species formulation allows the possibility 
X^s ds^R^s = 0, that is, the charge density can be zero in 
the rest frame of the fluid. 

In Eq.21, we have already identified the charge 
weighted flux sum (over matter species) with the cur- 
rent that drives the gauge field equation(14). This 
identification is perfectly sensible and is consistent with 
the choice of the fiuid equation of motion 

dt^T^ = gsnR.sU\.s. (22) 
Summing Eq.22 over species, we obtain 

^ d^T^ = a^r^"^ = dsFr'^r = F^J^.. (23) 

s s 



where T^" is the total fiuid tensor. Combining it with 
(I17|l . we arrive at the expected conservation law for the 
total energy momentum tensor ( matter plus field) 

9^(e^,+V)=0, (24) 

justifying the expression for the current. For the rest of 
the paper, we shall drop the species index unless it is 
essential for clarity. 

From the continuity equation, generalized to the non- 
Abelian case, 

P^F", = 0. (25) 

or 

d^{nRUn = -gnR,C''^'A\U^r (26) 

and from the the definition of T^^ (for a perfect fluid), 
we have 

d^T^^" = d^p + md^ifnRUa'Van 

= a> + muRUa^'d^ifUan - gnRmfUa^Ab^U'^'C''''' 
= ^''p + mnRUa^,[V^'fU'']a (27) 

The second term on the R.H.S of Eq.27 is related to Sa'^'" 
given by eauation ((T^ : 

Ua^Sa''" = Nd-^f + Uaf,[D^fU'']a (28) 

Combining eauations l27|l and 128|l . 

^,,T^'•' - d-'p + muRiUa^Sa'''' ~ Nd'^f) (29) 

where N stands for the dimension of the gauge group 
and we have applied V^^Ua^i — —N as & natural general- 
ization of U^Ufj_ = — 1 in the Abelian case. Substituting 
appropriately from equations H20I21() and H23|l . we find 

in 

dy - NmuRd-f = gnR[Fa^'' + -^a^^lC/a^ (30) 

In analogy with the Abelian case j^, we define the 
new unified non-Abelian tensor 

Mr = Fr + jsr, (31) 

representing the matter and gauge field (including their 
interaction). Then, equation H3UI) becomes 

d'^p - NmuRd^f = guRM^Ua^ (32) 

Defining the entropy a from Eqn. (|19|) . in analogy with 
a = ln[{p/K2){f)^exp{-^)], Eqs.30-32 may be 
combined to yield 

Td'a = gM^Ua^ (33) 

For a homentropic fluid (a relevant limit for the QGP), 
the equation of motion becomes even simpler: 

777 

gM^\U,a = [Fa^'' + -SanUa^.^O (34) 
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We have just shown the existence of a unified "minin- 
imally" coupled potential for hot non-Abelian fluids 

Tn 

Q''a-K + jfU", (35) 
with its corresponding field tensor 

= d,Q: - a.Q^^ + .gc^^gi'g^ (36) 

Through A/^'' and Q^, we have put the non-Abelian flow 
field and the non-Abelian gauge field on the same foot- 
ing. The unification opens up an opportunity to apply 
the powerful machinery of gauge theories to the unified 
gauge- fiow field; this formulation complements the pre- 
vious work on the subject [3.11^. 

The equation of motion given byj^ is the analogue 
of the two equations that Jackiw |T3 finds for the con- 
tinuum version of Wong's equations. However, unlike 
Jackiw's generalization, which seems to couple Yang- 
Mills fields with a fiuid of non-Abelian Charges, our pro- 
posal provides a natural non-linearity within the coupled 
system since S°'^i, contains an interaction with the Yang- 
Mills fields through the connection A"^. This natural 
non-linearity provides a mechanism to bypass the theo- 
rem forbidding the existence of solitonic configurations 
in a plasma (3- 

IV. TOPOLOGICAL INVARIANTS OF THE 
YANG- MILLS "MAGNETO-FLUID". 

In order to explore any interesting consequences of this 
formalism for Yang-Mills fluids, let us return first to the 
Abelian formalism and results on helicity conservation. 

The spatial components of the equation of motion in 
the Abelian case may be spelled out as: 

Uo{-S"' + F°'') + U,{-S'' + F^') = (37) 

q q 

Since i^"* is just the electric field, let us call the combined 
factor (J^S^^ + -F°'), the fiuid-generalized electric field 

. Let us make the same prescription for the magnetic 
components. Then, since is just the relativistic factor, 
7, and Ui — jU, eq. (pTzjl corresponds to: 

jE + jU X B ^0. (38) 
An immediate consequence is the condition 

E-B^O. (39) 

Since the product of Alf^^ — ^Sf^i, + F^^, with its dual 

Mf_,i, = \ep,vXpM^P, is proportional to E ■ B, Ea. l|5^ 
demands 

^M^-'M^, = 0. (40) 



Since the left hand side of this equation is a boundary 
term, we can assume the existence of a K^^, which for 
electrodynamics, is conserved 

^^K^' = 0. (41) 

It is well known that K^^ gives rise to helicity conserva- 
tion and if^ is identified with the fluid field equivalent of 
the Abelian Chern-Simons vector A^j^J^^"^ (see 13] ). The 
quantity ^M'^^Mfj.i, represents the topological winding 
number (charge) of solutions to the fiuid field equations 
of motion and the fact that it is zero supports the widely 
held view that an ordinary electron positron plasma, for 
example, does not support stable, self confining knot like 
solutions. This is upheld by a virial theorem due to 
Shafranov which states that a static configuration of a 
plasma in isolation is dissipative. Recently in jlTi] it has 
been proposed that this no go theorem is circumvented by 
introducing non linear interactions. We shall now show 
that the generalization to the YM plasma overides this 
limitation and can support stable knot like solutions. 

For the non-Abelian case, the spatial part of equa- 
tion l|M|) is some what more complicated: 

Uo^'ijSa^' + Fa°') + ji J Sa'' + Fan = (42) 

To manipulate the Eq. H42|) , we have to face the question 
of factoring the non-Abelian four- velocity t/^". At this 
stage there is no immediate compulsion for a factoring 
out of the generators (charges) of the gauge group. In- 
stead, one may assume that U"" ^Ua^ = tr groupU ^ = 
—N {N being the dimension of the gauge group) implies 
the existence of a full non-Abelian fiow and that the ve- 
locity four-vector is normalized in each fiow. In terms of 
the fluid-generalized electric and magnetic fields, equa- 
tion (03) becomes: 

J2lEa+^U- X Ba = (43) 

a 

Because of the trace over the group indices, Ea. H42|) 
implies that unlike in the Abelian case, the product 
M'^^M.aiiu 7^ 0. Consequently, the non Abehan gener- 
alization of the topological charge is not necessarily zero 
opening up the possibility of non trivial topological struc- 
tures being supported by the non-Abelian fiuid plasma. 

However, we may define a generalized Chern-Simons 
vector, as 

dpC>^ = \M-^MaXp. (44) 

In terms of the connection Q**, the generalized Chern- 
Simons vector C/i is given by: 

C^' = Q\[M''^''' - l:e''''^PC'''''QbxQcp] . (45) 
6 

To analyze this topological term, we look at the gauge 
properties of this potential Q"^. Under a gauge trans- 
formation Qg, 

A'^^QgA^n-' --id,ng)n;\ (46) 
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it is not difficult to see that: 

q; = QgQ^n-^ ~ ^{d,ng)n-\ (4?) 

This means that the non-Abelian fluid velocity vector, 
f/°ju transforms covariantly 

U;, = QgU^n-' (48) 

under a gauge transformation and can hence be 

strictly identified with a non-Abelian gauge connection. 
An immediate consequence for our analysis is that the 
that we have written above will indeed give us a gen- 
eralized Chern-Simons invariant associated with Q"^. In 
addition, the Yang-Mills connection v4°^ will provide us 
with the standard Chern-Simons invariant. 

The important point is that the minimal coupling pro- 
cedure introduced in fT^ and generalized here to the 
non-Abelian case, allows us to link these two Chern- 
Simons terms through the fluid velocity vector. Unlike 
the Abelian case (where the divergence of the generalized 
helicity four vector for the combined system was forced to 
be zero), we have two topological quantities in the non- 
Abelian case: one coming from the combined fluid-|-YM 
case and one from the YM case with the fluid velocity 
vector tying the two together. 

To understand the topological implications of our 
result, consider the transformation of the quantity 

I 



where , JF, 5 are the duals of M, F, S respectively and 
the integral is over spacetime. Thus in the non-Abelian 
magneto fluid we may associate this non zero Pontryan- 
gin index with a non-Abelian magneto fluid helicity im- 
plying the existence of stable self confining non dissipa- 
tive solutions. In fact, the non triviality of the Hopf 
invariant ensures that fiux lines can be knotted and soli- 
tonic configurations are inevitable. These can have a 
number of consequences, such as the existence of glue- 
balls, as knotted solitons having the non-Abelian helicity 
as a topological quantum number, which may survive in 
the quark gluon plasma in the interior of a heavy ion 
collision or in the early universe. 

V. CONCLUSION. 

We have given the foundations of a consistent theory of 
non-Abelian fiuid field system, in which the flow field and 



^ f C^cPx under gauge transformations. For this we 
use the wedge product notation AAF = e'^^'' AiFjkd^x to 
write 

^ = ^ / ^M^^a; =-^1 Tr{Q ^dQ-'^Q^Q^Q) 
= ^ / ^^(Q A M + A Q A Q) (49) 

Using the transformation properties established for Q we 
get the gauge transformed Ig as Ig = gij / Tr{Q, AMg + 
iQ' A Q' A Q') where Mg = flgAm-^ 

Thus we have 



^9-1=1^ j TridQgn-^ A dQgi}-^ A dUgQ-^] (50) 

which is the second Chern Class which describes the 
winding number of the manifold. This implies that the 
integral of the invariant trM^'^ A4^'^ will lead us to just 
the one appropriate Pontraygin invariant for the Yang- 
Mills gauge group that is used for the dynamics. From 
the minimal coupling prescription we have given above, 
we can write: 



(51) 



I 

the gauge field are "unified " in a single minimally cou- 
pled gauge flow field. We have shown that this gives rise 
to a quantity which is the fluid field generalization of the 
non-Abelian Chern Simons term, and shown that knot- 
ted fiuid-field non-Abelian solitons may exist. We can, 
by using standard techniques in pure Yang Mills theo- 
ries, find explicit forms of these and look for phenomeno- 
logical signatures in the context of QGP. The formalism 
is simple and unified and should lead to new and interest- 
ing phenomenon such as non-Abelian Alfven waves and 
other non-Abelian counterparts of magnetohydrodynam- 
ics which may lead to new signals for collective flow in 
the QGP. Such studies are under investigation. 
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